function dl = computeDescriptionLength(mix, nparsover2, ndata)

valid_inds = find(mix.priors > 0);  % find all the valid components
dl = sum(log(mix.priors(valid_inds)))*nparsover2 + (nparsover2+0.5)*length(valid_inds)*log(ndata);
